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Abstract. The new NASA Vision for Exploration, announced by 
President Bush in January 2004, proposes an ambitious program 
that plans to return astronauts to the moon by the 2018 time 
frame. A recent NASA study entitled “Affordable Fission Surface 
Power Study” recommended a 40 kWe, 900 K, NaK-cooled, 
Stirling convertors for 2020 launch. Use of two of the nominal 5 
kW convertors allows the system to be dynamically balanced. A 
group of four dual-convertor combinations that would yield 40 
kWe can be tested to validate the viability of Stirling technology 
for space fission surface power systems. 

The work described in this paper deals specifically with the 
control system for the 5 kW convertor described in the preceding 
paragraph. This control system is responsible for maintaining 
piston stroke to a setpoint in the presence of various disturbances 
including electrical load variations. Pulse starting of the Free 
Piston Stirling Engine (FPSE) convertor is also an inherent part 
of such a control system. Finally, the ability to throttle the engine 
to match the required output power is discussed in terms of 
setpoint control. Several novel ideas have been incorporated into 
the piston stroke control strategy that will engender a stable 
response to disturbances in the presence of midpoint drift while 
providing useful data regarding the position of both the power 
piston and displacer. 

Index Terms — Free piston Stirling engine, stroke control. 

INTODUCTION 

K. an kam [4] et al. have investigated the dynamic properties 
of the Free-Piston Stirling Engine (FPSE). The FPSE stroke 
amplitude is a function of electrical load. As the load 
resistance is decreased the stroke amplitude will be lowered to 
a new stable equilibrium point or stall if the resistance is too 
small to maintain oscillation. As the load resistance is 
increased the stroke amplitude will increase to a new 
equilibrium point or hit the mechanical limits. Modeling the 
dynamic behavior of the FPSE in response to load resistance 
changes is the focus of this work. 

The FPSE used in this work will be capable of producing 
400 volts RMS at 5 kilowatts with a 22 mm peak-to-peak 
stroke at a frequency of 85 hertz. The engine will be started by 
applying a short AC burst at 85 hertz. Once the engine stroke 
is sustained the stroke control circuit will maintain the stroke 
amplitude at the setpoint value. The setpoint range is between 
1 1 and 22 mm peak to peak. 


A master controller will monitor all subordinate controllers 
including: the stroke controller, the hot oil controller, and the 
cooling loop. The stroke controller is responsible for 
maintaining the FPSE convertor piston stroke at the desired 
setpoint. The hot oil controller is responsible for delivering the 
hot oil that is used to provide heat to the hot end of the FPSE 
convertor and keep it at the desired temperature. The cooling 
loop controller is responsible for keeping the cold end of the 
FPSE convertor at a desired temperature. In the event there is 
a problem with any of these subordinate controllers the master 
controller will implement a safe shutdown of the entire 
system. A touch screen will provide the user interface for the 
control system. Process data will be displayed and operator 
functions can be performed via the touch screen. For example, 
the stroke setpoint can be changed via the touch screen. 

The objective of this work was to understand the dynamic 
characteristics of a free piston Stirling engine and to establish 
an acceptable method of controlling the piston stroke of the 
engine. In order to experiment with various control strategies, 
without damage to the actual engine, a dynamic model of the 
engine was developed. Initially a linear model is used and 
later extended to a non-linear model. The non-linear model 
was then used to analyze different control strategies. Once a 
control strategy was selected the implementation of the 
control strategy is discussed. The dynamic model is shown in 
Figure 1. 


I. THE FPSE DYNAMIC MODEL 



Film re 1: Elements of the FPSE Model 

The convertor is mounted on a heavy steel plate that is 
sized so that the vibration amplitude of the convertor casing is 
limited to approximately 0.001” (25 |Jm). The engine casing 
can be considered as ground as shown in Figure 1. The 
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relative spring and damping parameters: K !2 , K 2 i, C [2 , and C 2 i 
represent the engine thermodynamics and are determined from 
the HFAST2 thermodynamics code. Due to the temperature 
difference between the expansion and compression spaces the 
two spring coefficients and the two damping coefficients are 
not equal to each other as would normally be the case in a 
passive dynamic system. The difference between the spring 
coefficients gives the dynamic system the ability to generate 
power. The hot end temperature is maintained at 650 °K while 
the cold end is 325 °K. 

The spring and damping coefficients Kn, K 22 , Cn, and C 22 
are determined from the STIRDYN dynamics code to achieve 
operation of the engine at the proper frequency (85 Hz) and 
displacer phase angle (70°). 

The damping coefficient connecting the piston to the 
alternator circuit, C 32 , has the units Volt-second/meter (E-s/m, 
where E is the alternator generated voltage) and the 
corresponding coefficient connecting the alternator circuit to 
the piston, C 23 , has the units Newton / ampere (F/A). 
Conservation of energy requires that the mechanical power 
delivered to the alternator by the piston (F*V) equals the 
mechanical power absorbed by the alternator (E*i), i.e. 

T7 r ■ F E 

-F-V = E-l, or, = — . (1) 

i v 

II. Model Equations 

The differential equations describing the dynamics are given 
in equation (2). 
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(4c) 


C = [o 0 1 0 0 o], (4d) 

The state variables were selected as follows: x | =q (charge), x 2 


• • • • 

m n x d + c n x d + c l2 x p + k u x d +k n x p =0, (2a) 


= q (current), x 3 = x p (piston displacement), x 4 = X (piston 


• • • • • 

m 22 x p + c 21 x d + c 22 x p + c 23 q+ k 2l x d + k 22 x p = 0, (2b) 

•• • • 

q+ c 32 x p + c 33 q+ k 22 q = v. (2c) 

If the temperature ratio of the engine was 1.0, then the k 12 
would be negative and equal to k 2! . If the displacer rod area 
were zero and the temperature ratio of the engine was 0, then 
both spring terms would be zero. From the analysis the load 
resistance is estimated to be 31.053 Q. Equation set (2) was 
used to develop a state space model (SSM) of the engine 
dynamics. The form of the SSM is given in equation (3). 

[x] = [A]x + [B]u and [ v] = [C]x (3) 


velocity), x 5 = x d (displacer displacement), x 6 = Xd (displacer 
velocity). Note that the output voltage is a function of piston 
velocity and that c 33 is the total resistance in series with the 
alternator. This includes the armature resistance and the load 
resistance. Also, c’ 33 (4.436 D) is the armature resistance. 
These coefficients for the linear model described in (4) are for 
the following operating conditions: head temperature (T H ) = 
650°K, cold end temperature (T c ) = 325°K, piston stroke (x P ) 
= 11 mm, and c 33 = 35.489 O. 

III. Alternator Equivalent Circuit 


A, B and C are given as: 



Figure 2: Electrical Equivalent Circuit 
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Figure 2 shows an equivalent electrical circuit for the FPSE 
alternator. The components are described as follows: Rc is the 
control resistance which is varied by the stroke controller 
output in order to maintain stroke setpoint. R L is the load 
resistance. The alternator inductance is m 33 The tuning 
capacitance is k33 (co-m 33 =1/ (co-k 33 ). The voltage 
proportionality constant is c 32 (V-s/m). Notice that c 33 — c’ 33 + 
R l ||Rc- It is the value of c 33 that determines the amplitude of 
the piston stroke. 


IV. The nonlinear model 

The nonlinearities are incorporated into the model using 
HFAST and STIRDYN to establish the relationships between 
the coefficients and the piston displacement. That is, Cjk = 
f(x p ), kj k = f(x p ), and m ik = f(x p ). Since a 6 th order model is 
used there are 36 coefficients to be established. After 
establishing all 36 coefficients for a number of piston stroke 
levels the nonlinear model was developed and programmed 
into Simulink® as shown in Figure 3. 



V. The open loop characteristics 

There are two possible FPSE modes of operation. The first 
is stable and oscillatory and the other mode is stable but 
damped. With all other coefficients constant, c 33 determines 
the mode of the dynamic system. If c 33 = 35.489 II the system 
is stable and oscillatory. The eigenvalues (open loop poles) 
are located in the complex frequency plane as follows: - 
210.69±j486.34, -70±j563.12, 0±j534.07. The root locus and 
response are shown in Figure 4. 


changes in resistance. The poles closest to the imaginary axis 
will dominate the response as poles with large real parts result 
in a faster response and their affect on total response may be 
small compared to those close to the imaginary axis. 

At some point reducing c 33 will result the inability for the 
FPSE to maintain oscillation due to the internal magnetic 
forces and the engine will stall. The engine will stall when c 33 
<15 12.. Figure 6 represents the response to a load change with 
c 33 changed from 40 to 14.3 O after ten seconds. In reality, c 33 
is changed by changing the value of Rl 
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If c 33 is decreased from 35.489 Q to 23.996 Q the stroke is 
reduced from 1 1 mm to 5 mm. The root locus for both values 
of c 33 are shown if Figure 5. 

The resistance change from 35.489 to 23.066 'Q has resulted 
in a significant change in the real parts of the complex 
conjugate pairs. This shows the sensitivity of these poles to 


x P = 1 1 mm , T h = 650°K 
T c = 325°K,c 33 = 35.489 'O 

Figure 4: Root Locus and Response with c 33 = 41.87 12 

Figures 5a and 5b show the root locus and pulse start 
response respectively for c 33 = 35.489 £2. Figure 5c and 5d 
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show the root locus and pulse start response respectively for 
C 33 = 23.966 D. 
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Figure 5a: Root Locus 


Figure 5b: Pulse Start Response 
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Figure 5d: Pulse Start Response 


Figure 6 shows the effect of decreasing the load resistance 
to a point where the FPSE stalls. In this case a resistance value 
of 15 'O will cause the FPSE to stall. 
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Figure 6: Response to Load Change 


In order to establish the open loop dynamic characteristics 
of this engine a series of step response tests were conducted 
and a transfer function was established for each. Figure 7 
shows the results of the step response tests. These tests were 
performed by changing R L . Equation (5) shows the transfer 
function for four such tests. 


scheme. The output of the controller can change the 
total resistance via proportional pulse width 
modulation (PWM) switching of the control load. 
This keeps the stroke at the setpoint value. The duty 
cycle of the PWM signal will determine the effective 
value of R c as shown in Figure 2. 

The setpoint controller is an algorithm that adjusts 
the stroke setpoint to match the user load 
requirements. Therefore, a non-dissipative type 
control scheme is realized which essentially throttles 
the FPSE. A block diagram of the control scheme is 
shown in Figure 8. Control load current, user load 
current, and terminal voltage serve as inputs to the 
setpoint control algorithm. The output of the setpoint 
control algorithm is the stroke setpoint. The 
controller is capable of holding the piston stroke at the 
setpoint in the presence of both thermal and electrical load 
disturbances. Flowever, the controller cannot cause the piston 
stroke to increase if the user load resistance Rl, which is in 
parallel with the control resistance Rc. is so small. 

The controller also monitors the piston amplitude 


Response to Resistance Step Change 
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Figure 7: Load Step Response 

Notice in (5) that the pole location and the gain change with 
each step change. This is due to the nonlinearity of the system. 
Flowever, each response can be approximated with first order 
dynamics. This will be useful when tuning the controller as 
discussed in the next section. For example, if we constrain the 
stroke range between 5 and 1 1 mm we can use the first order 
dynamics described by the second transfer function in 
equation set (5). 

VI. THE CONTROL STRATEGY 

One method of holding the piston stroke at setpoint is to 
change c 33 using a closed loop control system. A displacement 
transducer is used to measure the piston stroke and serve as 
the feedback, or process variable, in the closed loop control 



Figure 8: Control System Block Diagram 


midpoint to ensure that if the midpoint drifts off center the 
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controller will reduce the setpoint to prevent the piston from 
hitting the mechanical limits. This is only possible if there is a 
piston stroke measurement in place. If voltage is used to 
control the stroke then this drift information is not present and 
over stroking could occur if the piston drifts off center. Since 
the output voltage is proportional to the piston velocity, the 
power output is reduced as the stroke amplitude is decreased. 
The controller should be able to respond to user load changes 
without significant overshoot which might result in the 
mechanical limits being reached. 


Then at 20 seconds the load resiatance was changed from 40 
to 50 D. Finally, after 30 seconds the load resistance was 
changed from 50 to 100 XI The response plots are shown in 
Figure 9. Notice the response when the load resistance is 
changed from 40 to 50 12 s. 

The controller holds the stroke disturbance to about 9% 
when the load resistance is changed from 40 to 50 12. 
Increasing the integral gain (move the controler zero to the 
left) will further minimize the disturbance amplitude. 


From equation set (5) it is clear that the dynamics are not 
integrating processes. Therefore, integral action must be 
included in the control law to ensure that there is no offset 
error in the final stroke amplitude. Integral action places a 
pole at the origin and a zero on the real axis in the left plane. 
The placement of the zero determines the amount of integral. 
Placing this zero is part of the controller tuning. 

The stroke setpoint is manipulated based on the user load 
demand. That is, as the user load is increased, load resistance 
is lowered, and the stroke setpoint is increased. As the user 
load is decreased the stroke setpoint is decreased. The result is 
a non-dissipative control which throttles the FPSE to match 
the user load. The control load current is minimized to 
conserve fuel. 


VII. CONTROLLER TUNING 

A proportional-integral (PI) control law was incorporated in 
the controller. The PI control law transfer function is given by 
equation (6). 

U(s ) _ K p( S + /V _ 8.625(5 + 7.69) 

E(s ) 5 s 

With the control law in place a series of closed loop load step 
changes were made. The FPSE was pulse started and after 10 
seconds the load resistance was changed from 100 to 40 X2. 



time (sec) 


Figure 10 shows how the controller modulates the control 
resistance Rc to hold the piston stroke to a specific setpoint in 
response to a load disturbance. In Figure 10, Rl was changed 
from 40 to 50 12. 
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Figure 10: Controller Response 

Figure 1 1 shows the response when a step change, from 
0.0075 to 0.01 1 meters, is made in the stroke setpoint ten 
seconds after a pulse start. 
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Figure 9: Load Disturbance Response 


Figure 11: Setpoint Response 
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VIII. CONCLUSIONS 

The dynamics of a FPSE can be modeled by a 6 th order 
nonlinear model. It is possible to maintain stable operation of 
a FPSE, at a desired stroke amplitude setpoint, by modulating 
a control resistance in parallel with a user load. Since the 6 th 
order model is a type 0 dynamic system, integral action is 
required to maintain an offset error of zero in the steady state. 

From Figure 5 it is clear that the least dominant pair of 
complex poles are the most sensitive to load resistance. A 
43% reduction in resistance engendered a 61% change in the 
real part of the least dominant poles but, only a 1 .4% change 
in the most dominant pair. 

The controller can be tuned based on a first order model of 
the transfer function between the stroke amplitude and the 
total resistance. 

The work presented in this paper was based on a linear 
model of the FPSE that is currently being designed and built 
by Foster Miller, Inc. Once the engine is built the model and 
these results will be compared to the system. There are 
certainly a number of opportunities for nonlinearities to 
produce results that may differ form those obtained with this 
model. The nonlinearities can be identified and most of these 
can be incorporated into the model. 
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